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Abstract 

The radiation at collision of high-energy particles is formed over a rather 
long distances and therefore is sensitive to an environment. In particular 
the smallness of the transverse dimensions of the colliding beams leads to 
suppression of bremsstrahlung cross section for soft photons. This beam- 
size effect was discovered and investigated at INP, Novosibirsk around 1980. 
At that time an incomplete expression for the bremsstrahlung spectrum was 
calculated and used because a subtraction associated with the extraction of 
pure fluctuation process was not performed. Here this procedure is done. 
The complete expression for the spectral-angular distribution of incoherent 
bremsstrahlung probability is obtained. The case of Gaussian colliding 
beams is investigated in details. In the case of flat beams the expressions 
for the bremsstrahlung spectrum are simplified essentially. Comparison of 
theory with VEPP4 and HERA data is performed. Possible application of 
the effect to linear e^e~ collider tuning is discussed. 
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1 Introduction 



The formation of the bremsstrahlung process of high-energy particles occurs with 
extremely small momentum transfers. In the space-time picture this means that 
the process takes place at the large (macroscopic) distances. The longitudinal 
length (with respect to the direction of the initial momentum) of formation of the 
radiation usually is called the coherence (formation) length If. For emission of a 
photon with energy u the coherence length is //(cu) ^ e{6 — uj)/'m'^u), where e and 
m is the energy and mass of the emitting particle ( here the system ^ = c = 1 is 
used). If the particle experiences some action in this length, the radiation pattern 
changes (in the case when the action is the multiple scattering of the emitting 
particle one observes the famous Landau-Pomeranchuk effect [Q). 

A different situation exists in the bremsstrahlung process at the collision of 
electron and electron (positron) in colliding beams experiments. The point is 
that the external factors act differently on the radiating particle and on the recoil 
particle. For the radiating particle the criterion of influence of external factors 
is the same both at an electron scattering from a nucleus and at a collision of 
particles. For the recoil particle the effect turns out to be enhanced by the factor 
e^/m^, which is due to the fact that the main contribution to the bremsstrahlung 
cross section give emitted by the recoil particle virtual photons with very low 
energy 



qo ~ -I r, (1-1) 



so that the formation length of virtual photon is 

^^M = ^/(^o) = (1-2) 

This means that the effect for the recoil particles appears much earlier than for 
the radiating particles. For example, the Landau-Pomeranchuk effect distorted 
the whole bremsstrahlung spectrum in a TeV range (for heavy elements) while it 
turns out that the action on the recoil particle can be important for contemporary 
colliding beam facilities in GeV range p[. 

There are a few factors which could act on the recoil electron. One of them 
is the presence of an external magnetic field in the region of collision of particles 
If the formation length of virtual photon L„ turns out to be larger than 
the formation length Ini^^) of a photon with energy u; in a magnetic field H than 
the magnetic field will limit the region of minimal momentum transfers, which 
will lead to a decrease of the bremsstrahlung cross section and a change of its 
spectrum. Another effect can appear due to the smallness of the linear interval / 
where the collision occurs in comparison with Ly{uj) (see ( pr2|) ). This was pointed 
out in [^. 

In the experiment P| devoted to study of the bremsstrahlung spectrum da^{uj) 
carried out in the electron-positron colliding beam facility VEPP-4 of Institute of 
Nuclear Physics at an energy e = 1.84 GeV, a deviation of the bremsstrahlung 
spectrum from the standard QED spectrum was observed. This was attributed to 
the smallness of the transverse size of the colliding beams. Theoretically the prob- 
lem of finite transverse dimensions was investigated in |^ were the bremsstrahlung 
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spectrum at e"*"e~ collision have was calculated to within the power accuracy (the 
neglected terms are of the order I/7 = m/e) . Later the problem was analyzed 
in [§, IP, [|l^ where the found bremsstrahlung spectra coincide with obtained in 

0- 

It should be noted that in [0 (as well as in all other papers mentioned above) 
an incomplete expression for the bremsstrahlung spectrum was calculated. One 
has to perform the subtraction associated with the extraction of pure fluctuation 
process. Let us discuss this item in some details. The momentum transfer q at 
collision is important for the radiation process (the cross section contains factor 

at q^ <^ m^). At the beam collision the momentum transfer may arise due 
to interaction of the emitting particle with the opposite beam as a whole (due 
to coherent interaction with averaged field of the beam) and due to interaction 
with an individual particle of the opposite beam. Here we are considering the 
incoherent process only (connected with the incoherent fluctuation of density) 
and so we have to subtract the coherent contribution. The expression for the 
bremsstrahlung spectrum found in contains the mean value < q^ >, while the 
coherent contribution contains < q >^ and this term has to be subtracted. We 
encountered with an analogous problem in analysis of incoherent processes in the 



oriented crystals [TT| where it was pointed out (see p. 407) that the subtraction has 
to be done in the spectrum calculated in 0. Without the subtraction the results 
for the incoherent processes in oriented crystals would be qualitatively erroneous. 

In Sec. 2 a qualitative analysis of the incoherent radiation process is given. 
In Sec. 3 the general formulas for the spectral- angular distributions of incoher- 
ent bremsstrahlung are derived. The incoherent bremsstrahlung spectrum for the 
Gaussian beams is calculated in Sec. 4 in the form of double integrals. In specific 
case of narrow beams (the size of beam is much smaller than the characteristic 
impact parameter) the formulas are simplified essentially (Sec. 5). The experi- 
mental studies of effect were performed with fiat beams (the beam vertical size is 
much smaller than horizontal one). This specific case is analyzed in Sec. 6, while 
comparison with data is given in Sec. 7. In Sec. 8 the possible application to the 
linear e"'"e~ collider tuning is discussed. 



2 General analysis of probability of incoherent 
radiation 

In this section we discuss in detail the conditions under which we consider the 
incoherent radiation. One can calculate the photon emission probability in the 
target rest frame, since the entering combinations u/e and 7?? (7 is the Lorentz 
factor 7 = e/m, '(9 is the angle of photon emission) are invariant (within a rela- 
tivistic accuracy) and a transfer to any frame is elementary. We use the operator 



quasiclassical method |12], [13]. Within this method the photon formation length 
(time) is 

e' e' Ifo 

If = = — 

ekv euj{l — nv) ( ' 
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_ 1 _ 2ee' _ ie'-fcSr 
qmin ojm'^ um^ 
C=l + ^^d\ e' = e-uj, (2.1) 

where = ev^ {v^ = (1, v)) is the 4-momentum of radiating particle, 7c = ec/rric, 
Ec is the energy of target particle in the laboratory frame, iric is its mass, is the 
energy of radiating particle in the laboratory frame; = [u, oon)) is the photon 
4-momentum, i) is the angle between vectors n and v. 

In the case when the transverse dimension of beam a is a 3> //o the impact 
parameters g < Qmax = ^/o contribute. One can put that the particle density in 
the target beam is a constant, so that the standard QED formulas are valid. Note 
that the value Qmax is the relativistic invariant which is defined by the minimal 
value of square of the invariant mass of the intermediate photon In the case 
when the characteristic size of beams is smaller the value Qmax the lower value of 
is defined by this size. 

In the target rest frame the scattering length of emitting particle is of order 
of the impact parameter g. This length is much smaller than the longitudinal 
dimension of the target '-/J (/ is the length of target beam in the laboratory 
frame). So one can neglect a variation of configuration of the beam during the 
scattering time. A possible variation of particle configuration in the beam during 
a long time one can take into account in the adiabatic approximation. 

Another limitation is connected with the influence of value of transverse mo- 
mentum arising from the electromagnetic field E = |E| of colliding (target) beam 
on the photon formation length. This value should be smaller than the character- 
istic transverse momentum transfer ^m^/( in the photon emission process: 

eElf aNc 1 'ie''~^c^r 



{a^ + cry)l^c my/C uCm? 
2aN, 1 2€'€r AaN.jrXle' 2) 



here a = 1/137, N^. is the number of particles in the target beam, cr^ and ay are 
the vertical and horizontal transverse dimensions of target beam. Note that the 
ratio 7// is the relativistic invariant. This condition can be presented in invariant 
form 

-% « 1' (2-3) 

2 

7 UJ TfL 

where x = — + v x HI, -u = — , Eq = — = 1.32 ■ lO^^V/cm. Since the main 

Eq e' e 

contribution to the spectral probability of radiation gives angles -(9 ~ I/7 (C ~ 1) 
this condition takes the form xju <^ 1. For the case xju'^l the condition (|2.3|) 
can be satisfied for the large photon emission angles C — 7^^?^ > {x/'^Y^^ ^ 1- 
Under these conditions the formation length If = l/o/C decreases as {x/'^Y^^- The 
same inhibition factor acquires the bremsstrahlung probability fl^ . 



We consider now the spectral distribution of radiation probability in the case 
X ^ 1 (this condition is fulfilled in all existing installations), so 

X ~ «Ar,7 ^} « 1. (2.4) 
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Only the soft photons {uj < ^) contribute to the coherent radiation ("beam- 

strahlung" ) while the hard photon region 3> is suppressed exponentially as 
it is known from the classical radiation theory. As it was mentioned in the soft 
photon region (cj < ^ e), the spectral probability of bremsstrahlung is sup- 
pressed by the factor (cu/ex)^/^ only. On the contrary, the spectral probability of 
the bremsstrahlung is negligible comparing with the beamstrahlung taking into 
consideration that the mean square of multiple scattering angle during all time of 
beam collisions is small comparing with the value 1/7^: 

2 02 < > Sa^N^Xl ^ ^ ,^ , 

72 < ^2 Hs _ c_j^j^ ^ 2.5 

where L is the characteristic logarithm of scattering problem (in the typical ex- 
perimental condition L ~ 10). 

It was supposed in the above estimations of beamstrahlung probability that 
the radiation formation length is shorter than the target beam length 

I U \ Uj I 

Besides it was supposed that one can neglect a variation of the impact parameter 
Q and therefore of the transverse electric field Ex(^) during the beam collision. 
This is true when disruption parameter is enough small 

A = , \ , « 1, = z, y) (2.7) 

So, we consider the incoherent bremsstrahlung under following conditions: 

X < 1, - < 1, A, < 1. (2.8) 

u 

3 Spectral-angular distribution of the 
incoherent bremsstrahlung probability 

In this section we derive the basic expression for the incoherent bremsstrahlung 
probability at collision of two beams with bounded transverse dimensions. 

We consider first the photon emission at collision of an electron with one 
particle with the transverse coordinate x. We select the impact parameter qq = 
\qq\ which is small comparing with the typical transverse beam dimension a but 
which is large comparing with the electron Compton length Ac (Ac <^ f?o ^ 
a). In the interval of impact parameters g = |r_|_ — x| > qq, where r_|_ is the 
transverse coordinate of emitting electron, the probability of radiation summed 
over the momenta of final particle can be calculated using the classical trajectory 
of particle. Indeed,one can neglect by the value of commutators [[pxij^ij]! = Sij 
comparing with the value p±g in this interval {p±g > mgo ^ 1). In this case the 
expression for the probability has the form (see [l^, Eqs.(7.3) and (7.4)) 

dw = \M{Q)\'^Wrir^)d^r^d^k, (3.1) 
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where _ 

M(q) = — / R(t)exp(ik'x(t))dt, k' = -k. (3.2) 

Here w,.(r_L)(i^r_L is the probabihty to find the emitting particle with the impact 
parameter ^ = r_L — x in the interval di^g = d'^r±, R{t) = R{p{t)), kx{t) = 
ut — krit) (for details see [|r^, Sec. 7.1). Integrating by parts in the last equation 
and taking into account that |q±(^)| < I /go ^ we find 

X f°° , Kit) , ie , ^ d Rivi) 

M{q) ~ - — — / exp{ik'vt)— TrVrdt ^ - — 3.3 

In ^ Loo ^ ' dt k'v{t) 2TTy/LU ^^dp± k'v ^ ' 

where 



= P — n(npj ~ e\y — n 

m 



/oo Q poo ( it\ I 

ex^{ik'vt)ci{Q, t)dt = -— / exp — VUg"^ + t'^)dt 
-oo OQ J—oo \lf J 

= ^--ft'i(^) = 2aqminCKi{gqminO-^ (3-4) 
If g If g 

for the Coulomb potential, Ki{z) is the modified Bessel function (the Macdonald 
function), R{p±) has a form of matrix element for the free particles: 

R{p^) =vt'{A + i(tB)^,- a ~ ^-_^(e*u), 

B ~ ^(e* X (n - u)); u = ^, C = 1 + u^ fc'^; = q^U. (3-5) 
Zee m 

where the vector e describes the photon polarization and the spinors ips and ips' 
describe the polarization of the initial and final electrons correspondingly. 

In the interval of impact parameters g < Xc the expectation value of operator 
< q\M~^M\q > cannot be written in the form ( p.l|) since the entering operators 
become noncommutative inside the expectation value. However because of the 
condition Ac -C o" in this interval Wr{r±) ~ Wr-(x) + 0(— ) one can neglect effect 
of the inhomogeneous distribution. For the same reason in the calculation of 
correction to the probability of photon emission, which is defined as the difference 
of dw{a) and the probability of photon emission in a inhomogeneous medium, one 
can extend the integration interval into region g < go- 

In this paper we consider the incoherent bremsstrahlung which can be con- 
sidered as the photon emission due to fluctuations of the potential V connected 
with uncertainty of a particle position in the transverse to its momentum plane. 
Because of this we have to calculate the dispersion of the vector m(^) with respect 
to the transverse coordinate q: 

< rriimj > — < nii >< nij >= J mj(r_|_ — x)mj{r± — x)wc{x)d'^x 

— J ■mi{rj_—x.)wc{'x.)d'^x J mj(r_|_ — x)wc(x)(i^x, (3.6) 

where tOc(x) is the distribution function of target particles normalized to the unity. 
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As a result we obtain the following expression for the correction to the prob- 
ability of photon emission connected with the restricted transverse dimensions of 
colliding beams of charged particles: 



dwi 



a d?k 



Tij{e,pj_,s,s')Lij, Ti 



d R*{p±) 
dp±i k'v 



dp 



Q)d^Q^ (^J mj(g)wc{x — Q)d^Q^ Wr{x)d^x. (3.7) 



- / mi[Q)Wc[yi 



Averaging over the polarization of initial electrons and summing over the po- 
larization of final electrons we find 



ee' 



2eu, , 4(eu)2 ^ 

[eiUj + UiCj) H — — UiUj + -^i^ij 



C 



{3.t 



Note, that one can choose the real vector e since only the linear polarization could 
arise in the case of unpolarized electrons. 

After summation in ( p.8|) over the polarization of emitted photon we have 



Tij = ^ iv5^J - ^u^uA , V=^ + j, ( = 1+ 



(3.9) 



Finally, averaging the last expression over the azimuth angle of emitted photon 
we obtain 

h rr^,^. rrr.^_.. 4(C - 1) 



T 



2ee 



-um,, u{o = v 



(3.10) 



Substituting the expression obtained into ( p.7| ) we find the correction to the 
probability of photon emission connected with the restricted transverse dimensions 
of colliding beams of charged particles: 



dwi 



a e' duo 



U{OF{u;,C)dC, 



(3.11) 



where 



F(^)(cu, = ^1 Kl{7]g) {wri^ +q)- Wr{^)) w,i^)d^xd^g, 
F^'\cu,C) = ^l{^l K^ir]g)^w,i^~Q)d'g^ Wri^)d\, (3.12) 



here i] = QminC- 

Using the integral 



Kl{r]g)gdg 



2 



Kl{rig) ~ Ko{vg)K2{iig) 



(3.13) 
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and integrating by parts we obtain 



7]^ 



-2 



\ 2 



dg 



where 



Ki(rig)—Wc{'x.— Q)d'^g] WrMd^x 



(3.14) 



(3.15) 



In the general case the axes of colhding beams are displaced with respect each 
other in the transverse plane by the vector xq with coordinates Zq^Uq. In this case 
we have to consider 

^i;,(x)^t/;,(x + xo), F(i'2)(a;, C) ^ F^^'^H^, C, xq), ^q)^^q + ^o) (3.16) 

The first term in the expression for F{uj, Q in ( p.l4|) coincides with the function 
F{(jjX) defined in [0, Eq.(13). The second (subtraction) term in ( p.l4|) which 
is naturally arisen in this derivation was missed in Eq.(13), [0 as it was said 
above. The expression ( 3.11| ) is consistent with Eq.(21.6) in the book |jl3| (see 
also Eq.(2.2) in |]TT|]) where another physical problem was analyzed. It is the 
incoherent bremsstrahlung in the oriented crystals. 

Below we restrict ourselves to the case of unpolarized electrons and photons. 
Influence of bounded transverse size on the probability of process with polarized 
particles will be considered elsewhere. 



4 Gaussian beams 

For calculation of explicit expression for the bremsstrahlung cross section we have 
to specify the distributions of particles in the colliding beams. Here we consider 
the actual case of Gaussian beams. Using the Fourier transform we have 



w[yi] 



(27r)2 
Wr(q) = exp 



(i^gexp(— ■iqx)iy(q); 



, Wc[,^.) 



exp 



^ t 2 2 , 2 



y y> 



(4.1) 



where as above the index r relates to the radiating beam and the index c relates to 
the target beam, and Ay {a^ and ay) are the vertical and horizontal transverse 
dimensions of radiating (target) beam. Substituting ( |4.1| ) into Eq.( |3.15| ) we find 



*(^) = J d^qexp{-iqQ)exp 

1 



4S? 



Qy 



IT 



AD 



' y 



2(a^ + Ai) 



(4.2) 



Below we consider the general situation when the axes of colliding beams are 
displaced with respect each other in the transverse plane by the vector xq with 
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the coordinates Zq, Uq. This displacement has essential influence on the luminosity. 
For the processes for which the short distances are essential only (e.g. double 
bremsstrahlung p|) the probability of process is the product of the cross section 
and luminosity. The geometrical luminosity per bunch, not taking into account 
the disruption effects, is given by 

C = iV,iV,$(xo), (4.3) 

where as above N^. and N^. are the number of particles in the radiating and target 
beams correspondingly. We will use the same definition for our case. Then we 
have 

dw^ = <l>(xo)(i(T^, dai = $"^(xo)(iwi, (4.4) 

where dwi is defined in Eq. (|3.11| ). 

We calculate first the function F^^\uj,Q in Eq.( |3.12| ) for the case of coaxial 
beams when xq = 0. Passing on to the momentum representation with the help 
of formula ( [4.1D we find 

F^'\ujX) = / m,(q)^,(q)F2 (j^^ qdqd^, (4.5) 

where r] = qminC is introduced in (|3.12|) , 

F2 (l^) = 7 / - exp(-zqe))rf2p, 

+ 1 

F2{x) = ln(a; + Vl + x^) - 1, g™„ = m^cu /Ae^e' , (4.6) 

XV 1 + x^ 

here value qmin is defined in c.m.s. of colliding particles. The function F2(x) 
encounters in the radiation theory. To calculate the corresponding contribution 
into the radiation spectrum we have to substitute ( |4.5| ) into Eq.( ^.ll| ) and take 
the integrals. After substitution of variable in 



we obtain the integral over ( in Eq.( p.ll 



h - ^ + ^1 expi-s'OdC = f{s) 



where 



^(v- 8s2)Erfc(s) + 4e-^' + 2Ei(-s^), (4.8) 
2s 

s = wrqmin, r'^ = cos^ ip + E'^ sin^ ip (4.9) 
Making use of Eq.(14) we find for the spectrum 

nh^i^y Jo ^ cos^ V9 + ^ sm Jo 

= -^^2 2 ! ^ 2 • 2 • (4-10) 

Qniin COS^ Lf + T^y^ Sill if 
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This formula is quite convenient for the numerical calculations. 

In the case xq 7^ we will use straightforwardly Eqs.( |3.1^ ) and ( |4.4| ). Taking 
into account (14.21) we have for the difference 



A«(xo) = ;^($-^(xo)F«(xo) - <l>-i(0)F«(0)) 



27r 



^ J Kl{7]g)exp [-glJ:l - fpjs^] {exp [-2f?,^oS^ - 2f?,yoSj] - Ijrf'f?, (4.11) 



rC 

where the function F'^^)(xo) is defined in Eqs. (|3.12|) , (|3.16| ). 
Using the Macdonald's formula (see e.g. [|I^, p. 53) 



roo 

2Kl{r]g) = / exp 
Jo 



2. V 



'r/2\ dt 



2t t 



(4.12) 



and taking the Gaussian integrals over and gy we get 



A«(xo) 



1 r^^^pH^Mi 



e Jo 



exp 



2v4 

y 



rfdt 

~2r' 



(4.13) 

For the correction to the cross section (see Eqs.(14) and ( 4.10|) ) we have corre- 
spondingly 

rfa« = ^--[/(^n-) + ^^^H-,xo)], (4.14) 



where 



J«( 



u;,xoj 



f/(C)A«(xo)rfC. 



(4.15) 



Now we pass over to the calculation of the second (subtraction) term F^'^\uj, () 
in Eq.(|3A^). Using Eq.(|3) we get 

I = T] Ki{rig)-Wc{yi - Q)d^g 
J g 

where 



(27r)2 

S{q) = / Ki{r]g)— exp{iqQ)d'^ g 
J QP 



/ S{q)- exp{-iqx.)wc{q)(fq, (4.16) 

J (J 



21X1 / Ki(r]g)Ji(qg)gdg = 2Txi - 
J ri 



q 1 



rj q^ + rj^ 



Using the exponential parametrization 



1 



exp 



ds 



g2 _|_ ^2 4 _/q 

and taking the Gaussian integrals over and qy we obtain 



(4.17) 



(4.18) 



Jo 



exp 



r 



rfs 



X 



yey 

+ 2ay s + 2al, 



ds 



s + 2alJs + 2al 



(4.19) 



10 



where and ey are the unit vectors along axes z and y. Substituting ( [4.19|) into 
Eg. ( p. 121) , taking the Gaussian integrals over z and y and using Eq.( |4.4| ) we get 
the correction to the cross section 



(2) 



2a'^ e' duo 



(4.20) 



where 

j(2)(^,xo) = 

G(Si,S2,Xo) 



expizl^l + 2/oSj) ^ dsi ds2g I ^""'^^ j G(si, S2, xq 



1&2 



1/2 r 



V AB 



aia2 / 1 ZqU 



2+ A 



B 



'1 

2 5 



X exp 



■^(ai + a2)-f (61 + 62) 



Here the function g appears as a result of integration over C,: 



9{q) 



^ - 7 + :i ) exp(-g\^)^ = ( - ^ ) exp(-g 



f — 



4 



v-^] exp(-g^) - 2q^ 



( - -g^ ) Erfc(g) + -e"" + Ei(-g' 



In ( [4 .211 ) we introduced the following notations 
1,1 1 



2Ar 



2A2' 

y 



ai,2 



si,2 + 2a2- 



^1,2 + 2^2' 



A = ai + a2 + a, i? = 61 + 62 + 6, s = Si + S2. 



(4.21) 



(4.22) 



(4.23) 



5 Narrow beams 

This is the case when the ratio qmirj (S^ + T^y) ^ 1, so that the main contribution 
to the integral ( |4.1CI| ) gives the region g~(^~l, z ^ 1. Using the asymptotics 



of function F2{z) at 2; ^ 1 

F2{z) ~ \Yi{2zf - 1 (5.1) 

and the following integrals 
1 /-s^ dip 



271^2 Sj, Jo 2 cos2 + 2 sin^ ' 



2TrT,zT,y Jo Sj2 cQs2 _|_ ^-2 gj^2 ^ S^2 cQg2 _|_ 5]-2 gj^^a ^ 

/oo , , /-oo / 4 4 \ 

dq' (a - /5 In g^) ^ " ^ + ^ j exp{-q'C)dC 

^ - ^) [« + /?(2 + C)] + (5.2) 
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where C is Euler's constant C = 0.577..., we get for the function f^^^uj) (|4.14|) 
the following expression 



:r 2 In 



+ 3 + C +-, qmin<^{^z + ^y)- (5.3) 



This expression agrees with Eq.(24) of [0]. 

Under the assumption used in (|5.3|) and the additional condition qmin{zo + 
yo) -C 1 the main contribution to the integral in ( [4. 131 ) gives the region t 3> r/^. 
In this case one can use the asymptotic expansion Ki{z) ~ l/z{z -C 1). Then we 
have for the function J*^^)(ti;, xq) in Eq. ([4.14|) the following expression 



J 



CUjXoJ 



exp 



2 

; 

3 



t 



J 



+ 



dt 



t 



(5.4) 



The expression ( p.4| ) is consistent with Eq.(26) of [0. 

In the case (xq + o"^ + cry)qmin 1 ^^e main contribution to the integral in 
( |4.21| ) gives the interval ~ (xq 

of expansion over in Eq. (|4.22|) we get 



~ iXn + cr^ + (^y)qmin ^ 1- Keeping the main term 
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qmin\^ 



~ f . 

3 



(5.5) 



The same result can be obtained if one neglects the term containing rj'^ in the 
exponent of integrand in Eq.( |4.19| ). 

Summing the cross section da = dai^^ +da^i^ with the standard QED bremsstrahlung 
cross section 



da, 



2or e' duj 







V 



In 



m 

Qmin 



- 1 



we get the cross section for the case of interaction of narrow beams 



(5.6) 



da^ = daQ + dai 



2a-' e' duo \ 

m? e uj \ 

e e' 
v = - + - 



2 

V 

3 



21n- 



m 



+ C + 2 



e — UJ, 



(5.7) 



where J is given in ( |5.4| 



J_ = — — exp(2;oS^ + ?/oSj) / dsx \ ds2G{si, S2,xo), (5.8) 

l^z^y Jo JO 



where entering functions defined in Eqs.( [4.2ll) and ( [4.23| ). 

In the case of coaxial beams xq = 0, J = one can take the integral in (|5l 
over one of variables (for definiteness over S2) using the formula 



dx 
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bzby + hy/a^oTy 



(5.9) 



After this we have the simple integral over s = Si 



J40) = ^l+6,^l+6y{j,+ jy), 

i-oo I 

J^^y = / D^^y{s)dS, D^^y = J= ■ , (S.IO) 

where 

s 

a,,y = s{l + 6,,y) + 2aly{2 + 6,,y), 6.,, = + 1 + 5.,,, = (5.11) 

z,y z,y 

The cross section (5/7) differs from Eq.(24) of [0 because the subtraction term 
J_ is included. Without this term, generally speaking, the bremsstrahlung cross 
section would be qualitatively erroneous. In particular an appearance of the term 
J_ violates, generally speaking, the symmetry of radiation cross section in opposite 
directions in e~e~ {e~e~^) collisions. 

To elucidate the qualitative features of narrow beams bremsstrahlung process 
we consider the case of round beams where the calculation becomes more simple: 

(r, = ay = a, A, = Ay = A, = = = ^^^^^ ^2y 

2 

b = a, 6i,2 = ai,2, B = A, 5 = ^. (5.12) 

We consider first the case of coaxial beams (xq = 0, J = 0), 

ds 



J- = (1 + 6) 



[sil + 5)+2 + 6][s6 + l + 6] 



In the limiting cases the function J_ has the form 

J_((5»l)~^, J_((5 = l)=21n4 j_(5«i)^inl. (5.14) 
6 2d 

In the first case the subtraction term J_ is small. For the beams of the same size 
the subtraction term J_ contributes to the constant entering into the expression 
for the cross section. The subtraction term J_ modifies essentially the cross section 
in the case when the radius of target beam is much smaller than the radius of 
radiating beam. In this case the cross section ( [5.7| ) contains the combination 

— — — J_ ~ m m — ■ 

4E2 2 2(t2 

So in the all cases considered above the cross section defines the transverse di- 
mension of target beam. 

When the axes of round beams are displaced with respect each other in the 
transverse plane the integral in ( |5.4| ) is 



lib lib 2 

In — — - — J_ ~ m — m — - = myma) . (5.15) 



J 



dx 

Ei(d) - C - In ci. 



^ x + l 
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It is convenient in this case to calculate the function J_ using straightforwardly 
Eg. ( [4. 191) where we omit the term with rf in the exponent of integrand 



I. 



Q / exp 

/o 



ds 



2(t2 



2a 



2\2 



exp 



(5.17) 



Substituting this expression (I is defined in Eg. ( [4. 161) ) into the subtraction term 
Eg.( |3.12D and using the exponential parametrization 



1 r'^ 

— = / exp(— ^^s)(is 
Jo 



we obtain 



J_ 



ae 



ae 



d—di 



J ds J d'^ Qexi){—Q^s) 
dia \ 



1 — exp 



2^2 







1 



s + a + a~ 



s + a 
exp ( di 



exp 



+ 
a 



s + a + (7-72 



exp(-a(^ + xo) 



exp di 



s + a + a-V2^ 



+Ei di 



a2 + 2A2 



s + a + cr~2 
di = oXq 



ds 



ae 



Ei{di) - 2Ei di 



4 + vl 

2A2 ■ 



a2 + A2 ; 
(5.18) 



In the limit c/i — ^ the last expression goes over to Eg.( |5.13| ). 

When the displacement of the axes of colliding beams is large enough (xq ^ 
0"^ + A^) one use the asymptotic expansion of the function Ei(2;) in ( |5.18| ): 



Ei(2) ~ — (^1 + ^ ) , z > 1. 



(5.19) 



In this case the main terms in the difference J — J_ in Eg. (|5.71 ) are canceled: 

a 1 \ 2e'^a2 



= . 

d \d d\) (i Xq 



(5.20) 



The compensation of the main terms in ( p.l9| ) is due to the fact that the incoherent 
scattering originates on the fluctuations of the potential of the target (scattering) 
beam. Correspondingly we have for the mean sguare of the momentum transfer 
dispersion at the large distance from the target beam 



< q^(e) > - < q(e) >^oc 



1 



^ (xo + Qf xg 



'4(xoe) 



(5.21) 



Substituting (|5.20|) into Eg. (|5^ ) and multiplying the result by the luminosity 
(El): 

C = N^Nr— exp(-xgS=^) (5.22) 



IT 
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we have for the probabihty of bremsstrahlung of round beams moving apart at 
large distance 



a 



71 e u 



e' duj 



X 



exp(-XoS )ln — 



V - 



2^2 



2 
3 

+ 0(exp(-x2s2)) 



2(A2 



^0+^0 »1 
2(A2 + (t2) ' 



9L„(^o+l/o)«l-(5.23) 



According to ( [5. 231 ) when Xq increases so that one can neglect the first term in 
square brackets, the probability of bremsstrahlung of the round beams diminishes 



as a power of distance between beams (oc a /xq). The cross section Eq.( ^77|) in 
this case grows exponentially as e'^/d^. Let us note that without the subtraction 
term one has erroneous qualitative behaviour of probability (oc 1/xq). These 
circumstances explain also Eq. ( p-lSp for the coaxial beams: at integration over 
d'^Q the region contributes where < q^(^) > — < q(^) >^oc 1/ Q^, so that Q<a. 

Let us consider now the general case ^ for enough large displacement 
of beams Xq ^ ^- ^^lis case the main contribution into the integral I(x) 
(for rf = 0) in Eqs.( |4.l6| ),( t4.19D at large |x| ~ |xo| (see Eq.( |3.12| )) are given by 
large values s ~ Xq a'^y. Expanding the integrand over the powers cr'i^y/s and 
keeping after integration the two main terms of the decomposition over 1/x^ we 
get 



1 + 



r2\2 



(5.24) 



Expanding the function l/(xo + over the powers ^/xq at the integration over 
^ = X — Xq in Eq. (|3.12D ) we find 



/i2(xo + o^,(i)rf2e^l 

J Xn 



1 + 



+ 



(5.25) 



In this case the region t ~ 1/xq <^ contributes into the integral J Eq.(5.4)) 



Expanding the integrand over the powers tS^ ^ and keeping the two main terms 
of decomposition over 1/xq we have 



J 



^0 



2 2,2 



(5.26) 



For the difference J — we obtain finally 

1 



J - J_ 



exp(zoS^ + yoSj^ 



(xg) 



2\2' 



(5.27) 
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6 Narrow flat beams (cr^ ^ cr^, ^ 

Let us begin with the coaxial beams. We consider first the case where the size 
of radiating beam is much larger than size of target beam {5z^y <^ 1). In this 
case one can neglect the terms oc 5z,y, cr"^, in the functions a^^y and bz^y in the 
integral in Eq. (|5.10| ). Within this accuracy 



a,~s + 4aj, 6,-^ + 1, 6,-1. (6.1) 



After substitution in the integral Jy in Eq.( |5.10| ) s 4cr^s one gets 



ds (yf, 



Jyi^^) = / , ^ =-, K=^. (6.2) 



After substitution in the integral in Eq.( |5.10 ) s — 2A^/s one g so 
that 



J_(k) = 2Jl + SzJl + SyJy{K) ~ 2Jy{K) 



J_(/t < 1) ~ In-, J_(fi; > 1) ~ vrW-. (6.3) 

K V 

It is seen from the last equation that at A^ -C Cy the contribution of the term J„ 
into the cross section Eq. (|5.71) is relatively small. In the opposite case A^ ^ ay 
this contribution leads to change of the logarithm argument in Eq.( |5.7| ) 

2 In - In - ^ 2(ln(v^mA,) - ln(2y2^)) = 2 In (6.4) 

(L^ + Ly) K ay 2 

This is a new qualitative result. 

In the opposite case when the size of radiating beam is smaller or is of the 
order of size of target beam {5z^y > 1) the contribution into the integral Jz in 



Eq.( |5.10|) gives the region s ~ and into the integral Jy the region s ~ a. 



2 

y 

y 



Performing in the integral Jz the substitution s '^cr'^s and in the integral Jy the 
substitution s 2cr^/s one gets 



J. ^ 



°° ds 



2 + Syay Jo {{s + l)6z + l)Js{l + 6z) + 2 + 6. 



2 Az 1 



arctan 



7 ~ 

a. 



2 + SyCry ^6,{2 + 6z] 

A, r°° ds 



yJo ^^S + l)i6y+l) + S)^iS + l)6y + 



2 Az 1 
— arctan ■ 



2 + Sy(^y J6y{2 + 6, 



j_ = ^i + s^^i + Sy(j^ + Jy) 

^ ?^i^£±5^ Letan ^ ' + arctan ^ ' ] .(6.5) 

^2TTy ^y V ^Sz{2 + 6z) ^5y{2 + 5y)^ 
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In the case S^^y ^1, <^ this formula is consistent with Eq.( |6.3|) . 

Now we go over to the case of the displaced beams. For enough large displace- 
ment of the beams the formulas ( ^.71) and (|5.27|) are valid. So the intermediate case 
is of interest. As an example we consider the case ^ ^ + A^, Vq CTy. 
In this case the contribution to the integral in (|5.4| ) gives the interval <^ t ~ 
Zq"^ <^ S^. Keeping the main terms of decomposition over 

tS;2 ^ I and > 1 

we have 

J^^r exp {^El - 4t) 4 = ^ exp(.o^S^). (6.6) 

Under these conditions (xq ^ cr^) the contribution into the integral for J_ in 
(|5.8| )of the term in the function G(si,S2,xo) Eq.( |4.21| ) in the square brackets 
containing 6162/-S is defined by the function Jy in Eq.( p.5|) to within the terms 
~ Zo/ay. In the term containing aia2/A (which we denote by J^^) the main 
contribution gives the summand ZqO^/A^ in the interval cr^ ^ Si,2 ~ ^0 ^ '^2 
where 

11 11 
si,2 2(T^ 2A^ 



As a result we obtain 



A, VI + 5. a/1 + ^ 

TT 



y 



l2 + 5y 



J-J--Jje''^h{zo), d, = zlJ:^ 



Z1 



Hzo) = 1 - V^EiM^ 1+2 Cretan , ' 1 . (6.8) 

Jm + 6y) V ^ J6y{2 + 6y) 



It should be noted that for the flat beams the probability of radiation as a function 
of distance between beams (for the considered interval) decreases more slowly 
oc \l\fdz than for the round beams given in Eq.( |5.23|) 



dwi^^^N^N:^\lY.Jly'--{v-\ 



(6.9) 



Compensation in the difference J — J- begins in the region zq ^ ay + Ay were 
Eq.(|6.8D is not valid and one have to use more accurate Eq.(p.8|). In the region 
zq ^ ay + /\y the probability of radiation decreases as l/zg according to Eqs.( [4.4|) , 
( |5.7| ), ( |5.27| ) provided that one can neglect the exponential term in the square 



brackets in Eq. (|6.9|) (compare with Eq.( ^.23|) ) 



dw!;{zo)^2N,N,— '-^-(v-^]—, zo:^yo. (6.10) 
' Tc Zq e \ 3J u 



^Xlale' { 2\ duj 
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7 Observation of beam-size effect 



Above we calculated the incoherent bremsstrahlung spectrum at collision of elec- 
tron and positron beams with finite transverse dimensions. This spectrum dif- 
fers from spectrum found previously in 0, P] because here (in contrast 
to previous papers) we subtract the coherent contribution. In general expres- 
sion for correction to the probability of photon emission (|3.11|) the subtraction 
term is F^'^\uj,Q- For the coaxial beams for numerical calculation it is conve- 
nient to use Eqs. (14.101 ), (|4.2CI| ) and (|4.21| ). In the last equation one have to put 
Vo = Zo = 0- In the case of collision of narrow beams the subtraction term in 
the bremsstrahlung spectrum (|5.7| ) is J_. The dimensions of beams in the exper- 
iment 1^ were cr^ = = 24 fj,m, ay = Ay = 450 fim, so this is the case of flat 
beams. The estimate for this case (|0| ) gives J_ ~ {A/3^/3)7raz/cry <^ 1. This 
term is much smaller than other terms in ( [5.71 ). This means that for this case the 
correction to the spectrum calculated in ^ is very small. 

The result of calculation and VEPP4 (INP, Novosibirsk) data are presented in 
Fig.l where the bremsstrahlung intensity spectrum uda/du is given in units 2arQ 
versus the photon energy in units of initial electron energy {x = uj/e). The upper 
curve is the standard QED spectrum, the three close curves below are calculated 
using Eqs. ([4.10|) and ( [4.20| ) for the different vertical dimensions of colliding beams 
(equal for two colliding beams a = az = A^):^ = 20 fim (bottom), cr = 24 /xm 
(middle), cr = 27 fj,m (top) (this is just the 1-sigma dispersion for the beams 
used in the experiment). We want to emphasize that all the theoretical curves 
are calculated to within the relativistic accuracy (the discarded terms are of the 
order m/e). It is seen that the effect of the small transverse dimensions is very 
essential in soft part of spectrum [sX uj/e = 10~^ the spectral curve diminishes in 
25 %), while for uj/e > 10~^ the effect becomes negligible. The data measured in 
are presented as circles (experiment in 1980) and as triangles (experiment in 
1981) with 6 % systematic error as obtained in (while the statistical errors are 
negligible). This presentation is somewhat different from [^. It is seen that the 
data points are situated systematically below the theory curves but the difference 
is not exceed the 2-sigma level [Q. It should be noted that this is true also in the 
hard part of spectrum where the beam-size effect is very small. 

The last remark is connected with the radiative corrections (RC). The RC 
to the spectrum of double bremsstrahlung |T^ (this was the normalization pro- 
cess) are essential (of the order 10 %) and were taken into account. The RC to 
the bremsstrahlung spectrum [|T^ are very small (less than 0.4 %) and may be 
neglected. It should be noted that the RC to the bremsstrahlung spectrum are 
insensitive to the effect of small transverse dimensions. 

The dependence of bremsstrahlung spectrum on beams characteristics was 
measures speciflcally in ||^. The first is the dependence of bremsstrahlung spec- 
trum on vertical sizes of beams cx^. It is calculated using Eqs. (|4.10| ) and ( |4.20| ) 
for uj/e = 10"'^. The result is shown in units 2arQ in Fig. 2. The data is taken 
from Fig. 7 in 0. The second is the measurement of dependence of bremsstrahlung 
spectrum on the vertical displacement of beams zq. It is calculated using Eqs. (|5.4|) 
and ( ^.8| ) for uj/e = 10~^. Because of displacement it is necessary to normalize 
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the spectrum on the luminosity 



TX 

see Eq.( |4.3| ). This means that when we compare the bremsstrahlung process 
(where the beam-size effect is essential) with some other process like double 
bremsstrahlung used in (which is insensitive to the effect) we have to mul- 
tiply the cross section of the last process by the luminosity C This is seen in 
estimate Eq. (|6.9| ): after taking out the exponent e~'^'= we have the luminosity as 
the external factor and in expression for ratio N^/N2^ (which was observed in 0) 
the cross section of double bremsstrahlung will be multiplied by the luminosity. 
After this operation the second term in square brackets will contain the combina- 
tion e'^'' h{zQ) / \/d^ which grows exponentially with the displacement zq increase. 
The normalized bremsstrahlung spectrum is shown in units 2arQ in Fig. 3. So, 
the very fast (exponential) increase with zq is due to fast decrease with zq of the 
double bremsstrahlung probability for the displaced beams. The data is taken 
from Fig. 8 in [||. It should be noted that in soft part of spectrum the dependence 
on photon energy u is very weak. It is seen in these figures that there is quite 
reasonable agreement between theory and data just as in This means that 
contribution of J_ term which is calculated only in the present paper is relatively 
small. 

One more measurement of beam-size effect was performed at HERA electron- 
proton collider (DESY, Germany) [|18|. The electron beam energy was e=27.5 GeV, 
the proton beam energy was £p=820 GeV. The standard bremsstrahlung spectrum 
for this case is given by Eq. ([5.6|) where qmin should be substituted: 

D _ ^m^rrip 

Qmin ~^ Qmin ' ['■^) 

here is the proton mass. In this situation the characteristic length is /j^ = 
l/q^in the photon energy u; = 1 GeV one has ~ 2 mm. Since the 

beam sizes at HERA are much smaller than this characteristic length, the beam- 
size effect can be observed at HERA. The parameters of beam in this experiment 
were (in our notation): az = = (50 ^ 58)/im, (Xy = Ay = (250 -i- 290)//m. In 
part of runs the displaced beams were used with zq = 20 fim and yo = 100 fim. 
The bremsstrahlung intensity spectrum ujda/duj in units 2arQ versus the photon 
energy in the units of initial electron energy [x = uj/e) for the HERA experiment 
is given in Fig. 4. The upper curve is the standard QED spectrum. We calculated 
the spectrum with beam-size effect taken into account for three sets of beams 
parameters; the set 1: = = 50 jj,m, (Xy = Ay = 250 /im, zq = yo = 0, 
the set 2: az = A^ = 50 fim,ay = Ay = 250 jj,m,ZQ = 20 ^m,yo = 0, the 
set 3: az = Az = 54 fj,m, Oy = Ay = 250 /im, zq = y^ = 0. The result of 
calculation is seen as two close curves below, the top curve is for the set 3, while 
the bottom curve is actually two merged curves for the sets 1 and 2. Since the 
ratio of the vertical and the horizontal dimensions is not very small, the general 
formulas were used in calculation: for coaxial beams Eqs. ( ^.11| ) and ( |4.20| ),and for 



displaced beams Eqs. (|4.14|) and (|4.20| ). It should be noted that the contribution 
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of subtraction term (Eq.( [4.20D ) is quite essential (more than 10%) for the beam 
parameters used at HERA. The data are taken from Fig.5c in [Q. The errors are 



the recalculated overall systematic error given in |T8[. It is seen that there is a 
quite satisfactory agreement of theory and data. The final data are given in [|18| 
also as the averaged relative difference 6 = {daqED — da^s) I daqEo (where daqED 
is the standard QED spectrum, a^s is the result of calculation with the beam-size 
effect taken into account) over the whole interval of photon energies (2-8 GeV), 
e.g. for the set 1 5^^ = (3.28 ±0.7)%, for the set 2 4^ = (3.57±0.7)%, for the set 
3 6ex = (3.06 ± 0.7)%, [g. The averaged <5> over the interval 0.07 < a; < 0.28 
(or 1.95 GeV< u < 7.7 GeV) in our calculation are for the set 1 is < 5 >=2.69%, 
for the set 2 is < 5 >=2.65%, for the set 3 is < 5 >=2.54%. So, for these data 
there is also a satisfactory agreement of data and theory (at the 1-sigma level, 
except set 2 where the difference is slightly larger). 

So, the beam-size effect discovered at BINP (Novosibirsk) was confirmed at 
DESY (Germany). Of course, more accurate measurement is desirable to verify 
that we entirely understand this mechanism of deviation from standard QED. 



8 Conclusion 

Above the influence of the finite transverse size of the colliding beams on the inco- 
herent bremsstrahlung process is investigated. Previously (see papers 0, 0, 0, 



Top for analysis of this effect an incomplete expression for the bremsstrahlung 
intensity spectrum was used because in it the subtraction was not fulfilled. It 
is necessary to carry out this subtraction for the extraction of pure fluctuation 
process which is just the incoherent bremsstrahlung. We implement this pro- 
cedure in the present paper. We indicated the cases where the results without 
the subtraction term are qualitatively erroneous. The first this is the case when 
the transverse sizes of scattering beam are much smaller than the corresponding 
sizes of radiating beam. For coaxial round beams see e.g. Eq.( ^.15| ) and for flat 
beams Eq.( |6.4| ). In contrast to previous papers here we draw a conclusion that the 
bremsstrahlung cross section is determined by the transverse sizes of scattering 
beam. 

The new qualitative result is deduced for the case when the displacement 
of beams is enough large. Then the square of momentum transfer dispersion, 
which determines the bremsstrahlung cross section, decreases with displacement 
increase faster than mean value the momentum transfer squared (see Eqs.( ^.2lD , 
(|5.27| )). As it was noted in Sec. 7, it is necessary to normalize the spectrum on the 
luminosity for displaced beams. Then the bremsstrahlung cross section grows ex- 
ponentially with displacement increase. This very fast (exponential) increase with 
Zq is due to fast decrease with zq normalization process probability for displaced 
beams. 

For Gaussian beams the expression for the bremsstrahlung spectrum is ob- 
tained in the form of double integrals convenient for numerical calculations (see 
Eqs. (^4.10| ), ( [4.20|) and (|4.21|) ).For soft part of spectrum we deduced the general ex- 



pression for spectrum which is independent of minimal momentum transfer g. 
and is defined only by transverse size of beams (see Eqs.(p.3|), ( |5.4| ) and (jST 
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The important feature of the considered beam-size effect is smooth decrease of 
radiation probabihty with growth of displacement of beams. For the flat beams 
we see in Eqs. ( |6.9| ), ( |6.10| ) that the main (logarithmic) term in expression for the 
probability decreases exponentially (oc exp(— 2;qS^) as luminosity), but there is the 
specific long-range term (x I/zq which results in quite appreciable radiation prob- 
ability even in the case when beam the displacement is large. This phenomenon 
may be helpful for tuning of high-energy electron-positron colliders. As an exam- 
ple we consider the "typical" collider were the beam energy is e = 500 GeV and 
the beam dimensions are equal and = 5nm and ay = lOOnm. The beam-size 
effect in this collider is very strong and for x = 10~^ the intensity spectrum is only 
~ 0.3 part of the standard udaQEoi!^) / du. The dependence of bremsstrahlung 
probability on the displacement distance zq (in nm) is shown in Fig.5. It is cal- 
culated using Eqs. ( p.6| )- |578|) for soft photons with x = 10"^ (the asymptotic 
formulas ( |6.9| ) are ( |6.10| ) are not enough accurate in this case). Actually the de- 
pendence on photon energy is contained in the external factor (1 — x){v{x) — 2/3). 
The curve in Fig.5 reflects the main features mentioned above. One can see that 
even for zo=100 (2:0 = 20 a^) the cross section is ~ 0.002 part of very large 
bremsstrahlung probability at head-on collision of beams. So, measuring the radi- 
ation for displaced beams one can estimate magnitude of displacement of beams. 
This information may be useful for beam tuning. 
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Figure captions 



Fig.l The bremsstrahlung intensity spectrum uoda/du in units 2arQ versus 
tlie plioton energy in units of initial electron energy [x = oj/s) for VEPP4 
experiment. The upper curve is the standard QED spectrum, the three 
close curves below are calculated for the different vertical dimensions of 
colliding beams (equal for two colliding beams o" = cr^ = A^):^ = 20 /im 
(bottom), cr = 24 iim (middle), a = 27 /im (top). The data measured in 
are presented as circles (the experiment in 1980) and as triangles (the 
experiment in 1981) with 6 % systematic error as obtained in 

Fig. 2 The bremsstrahlung intensity spectrum uda/dw in units 2ar'^ versus 
the vertical sizes of beams cr^ (in /xm). The data taken from 

Fig. 3 The normalized to luminosity C the bremsstrahlung intensity spec- 
trum uda/dw in units 2arQ versus the vertical displacement of beams Zq (in 
/im) . The data taken from . 

Fig. 4 The bremsstrahlung intensity spectrum ujda/duj in units 2ar'^ versus 
the photon energy in units of initial electron energy (x = oj/e) for the 
HERA experiment. The upper curve is the standard QED spectrum, the 
two close curves below are calculated with the beam-size effect taken into 
account: the bottom curve is actually two merged curves for sets 1 and 2 
(the set 1 is cr^ = = 50 /im, ay = /S.y = 250 /im, zq = uq = 0, set 2 is 
0"^ = A^ = 50 nm,ay = Ay = 250 fim,ZQ = 20 /im, ?/o = 0); while the top 
curve is for set 3 {az = A^ = 54 /im, ay = Ay = 250 fim, zq = yo = 0) . The 
data taken from Fig. 5c in |[T8 |. 



Fig. 5 The spectral intensity probability udw^/duj normalized to one particle 
in the beam in units 2arQT,zT,y/'7i versus the vertical displacement of beams 
Zo (in nm). 
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